Energy Landscape, Anti-Plasticization and Polydispersity Induced 
Crossover of Heterogeneity in Supercooled Polydisperse Liquids 
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Polydispersity is found to have a significant effect on the potential energy landscape; the average 
inherent structure energy decreases with polydispersity. Increasing polydispersity at fixed volume 
fraction decreases the glass transition temperature and the fragility of glass formation analogous to 
the antiplasticization seen in some polymeric melts. An interesting temperature dependent crossover 
of heterogeneity with polydispersity is observed at low temperature due to the faster build-up of 
dynamic heterogeneity at lower polydispersity. 
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Polydispersity is ubiquitous in nature. It 
is present in clays, minerals, paint pigments, 
metal and ceramic powders, food preservatives 
and in simple homogeneous liquids. It is com- 
mon in synthetic colloids, which frequently ex- 
hibit considerable size polydispersity []| and is 
also found in industrially produced polymers, 
which contain macromolecules with a range of 
chain length. Polydispersity has significant ef- 
fects on both the structure and dynamics of 
the system. Experiments [2| and simulations 
[3, 0] on colloidal systems show that increas- 
ing polydispersity, at a constant volume frac- 
tion, lowers structural correlations, pressure, 
energy and viscosity. Polydisperse colloidal sys- 
tems are known to be excellent glass formers. 
William et al Q suggest that colloidal glass for- 
mation results from a small degree of particle 
polydispersity. Crystal nucleation in a polydis- 
perse colloid is suppressed due to the increase 
of the surface free energy Q . Studies by several 
groups [3| have shown that the glass becomes 
the equilibrium phase beyond a terminal value 
of polydispersity. 

Despite being natural glass formers, relation- 
ships between polydispersity, fragility, energy 
landscape and heterogeneous dynamics have 
not been adequately explored in these systems. 
Because these systems exist in the glassy phase 
over a wide range of polydispersity, they offer 
opportunity to test many of the theories and 
ideas developed in this area in recent years. 
We find that polydispersity introduces several 
unique features to the dynamics of these sys- 
tems not present in the binary systems usually 
employed to study dynamical features in super- 
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cooled liquids and glasses. 

In this work we particularly investigate how 
polydispersity influences the potential energy 
landscape, fragility and heterogeneous dynam- 
ics of polydisperse Lennard-Jones (LJ) systems 
in supercooled regime near the glass transition 
[a|. The polydispersity in size is introduced by 
random sampling from a Gaussian distribution 
of particle diameters, a. The standard devia- 
tion 6 of the distribution divided by its mean 
a gives a dimensionless parameter, the polydis- 
persity index 5 = = . The mass rrii of particle i 
is scaled by its diameter as m^ = m{^)^. Mi- 
cro canonical (NVE) ensemble MD simulations 
are carried out at a fixed volume fraction, (j) on 
a system of TV = 864 particles of mean diameter 
(7 = 1.0 and mean mass m = l.O for S = 0.10, 
0.15 and 0.20 at = 0.52 and S = 0.10 and 
0.20 at = 0.54. All quantities in this study 
are given in reduced units (length in units of a, 
temperature in units of jf- and time in units 

of r = (^^^^^^-)2). The LJ interaction parame- 
ter e is assumed to have the same value for all 
particle pairs. 

At large supercooling the system settles into 
glassy phase. We first analyze the system from 
the perspective of potential energy landscape 
(PEL), which has emerged as an important tool 
in the study of glass forming liquids d, ^Q, [ll| . 
FiglBJa) and (b) show the variation of the aver- 
age inherent structure energy ((e/5)) with tem- 
perature (T) at both the volume fractions stud- 
ied. The value of {ejs) remains fairly insen- 
sitive to the variation in T at high T before 
it starts to fall with T (around T ^ 1.0). It 
has been established earlier in the context of 
the binary mixtures [9|, [l^l that the start of 
fall in (e/s) coincides with the onset of non- 
exponential relaxation in the time correlation 
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FIG. 1: (a) and (b) Temperature dependence of 
the average inherent structure energy, (e/s). For 
Fig[T](a), (b) & (c), fiUed circles, stars and trian- 
gles are for 5* = 0.10, 5" = 0.15 & S = 0.20 at 
= 0.52 and filled diamonds and squares are for 
S = 0.10 & 5 = 0.20 at <;A = 0.54, respectively. 
(c) The stretched exponent /3 vs. T obtained by 
fitting KWW equation to self-intermediate scatter- 
ing function, Fs{kmax,t) where kmax ~ 7.0. The 
lines are guide to the eye. Comparison between 
(a)/(b) and (c) shows that the fall of (e/s) corre- 
sponds to the onset of non-exponential relaxation 
in F,{kmax,t). (d) 13 vs. T from F,(fc,t) for dif- 
ferent k values. Data shown for S — O.l(Sl) & 
S = 0.2(S2) at (^ = 0.52. S = 0.15 omitted for 
clarity. 
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FIG. 2: (a) Average radial distribution functions 
(rdf) for the parent liquid (solid line) and the inher- 
ent structure (dashed line) for S = 0.10 (black) and 
S = 0.20 (red) systems at T = 0.50 and cp = 0.52. 
(b) Relaxation time r from KWW fit to Fs{kmax,t) 
for 5" = 0.10 (filled circles) and S = 0.20 (filled tri- 
angles) at (/) = 0.52. [Inset: Critical temperature T^ 
for particles of different sizes a^ obtained from the 
MCT equation, D' ~ (T-T^y for S = 0.10 (open 
circles) and S = 0.20 (open triangles) at = 0.52.] 



functions of the system. Wc show in Fig [TJc) 
that this correlation continues to hold in poly- 
disperse systems. The fall of (e/s) with T is 
consistent with the Gaussian landscape model. 



The average inherent structure energy de- 
creases with polydispersity (Fig[Ha) and (b)), 
which indicates that the packing is more effi- 
cient at higher S. In Fig [5] we plot the inherent 
structure (IS) and the parent liquid radial dis- 
tribution functions (rdf). At S = 0.20 there is 
hardly any difference between the rdf of the par- 
ent liquid and the IS. The coordination number, 
Nc at S = 0.10 and S = 0.20 obtained from 
the IS rdf are 13.1 and 14.6, respectively. This 
shows that packing is more efRcient at higher 
S and one would expect a slowing down of dy- 
namics at higher S. Instead, wc find that similar 
to colloidal hard spheres polydisperse LJ sys- 
tems also show a speed up of relaxation with 
S. The presence of smaller particles at higher S 
provides some sort of lubrication [ij, (Tj] , which 
speeds up the dynamics of the whole system. A 



plot of the Mode Coupling Theory (MCT) [3 
critical temperature T^ for particles of different 
sizes cr* (inset of Fig [H^b)) shows that the T^ 
for the largest-sized particles in S=0.20 system 
is smaller than the smallest-sized particles in 
S=0.10 system. This tells us that not only the 
smaller particles in S" = 0.20 system but the 
whole system has a faster relaxation. The rate 
of growth of relaxation time upon lowering of T 
decreases with S (Fig[21[b)). Hence as the sys- 
tem is cooled, vitrification is expected to occur 
at a lower T for the system at higher S. This 
should lead to a lowering of the glass transition 
temperature with S. 

Fragility is a term being used to character- 
ize and quantify the non-Arrhenius transport 
behavior in glass-for ming liquids as they ap- 
proach glass transition 24| . To study the effect 
of polydispersity on fragility, wc plot the diffu- 
sion coefficients in an Angell-like fragility plot 
in Fig[3l The plot clearly shows that increasing 
polydispersity at fixed volume fraction reduces 
the fragility of the liquid so that the system is a 
stronger glass former at higher polydispersity. 
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FIG. 3: Angell-like fragility plot at different S for 
the two cj> studied. The thick lines are VFT fit to 
the diffusivity data, D = Doeap( j_^ ). The refer- 
ence temperature Tr is chosen such that D{Tr) = 
4.5 X 10~^. The VFT extrapolation is used to locate 
Tr. The plot shows that fragility decreases with S 
and that for a given S fragility increases with in- 
crease in (f). [Inset: Strength parameter m (where 
m = -jP- QM) obtained from VFT fit as a function 
of S atA = 0.52.1 



This effect is analogous to the antiplasticiza- 
tion that has been observed in polymer melts 
|16| . PEL analysis shows that the antiplasti- 
cized system has smaller barriers to overcome 
in order to explore the configuration space [l3| • 
In the rest of the paper we explore the correla- 
tions between fragility and non-exponential re- 
laxation/heterogeneous dynamics. 
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FIG. 4: The non-Gaussian parameter, a2{t) for 
5 = 0.10 (solid line),^ = 0.15 (dashed line) and 
S = 0.20 (dot-dashed line) at four different T de- 
picting the crossovers between different S. Data is 
shown for S = 0.52. 



if we look only at the /? values and not its T- 
dependence we find that at high T, stretching is 
anti-correlated with fragility whereas at low T, 
we get the reverse scenario where the stretch- 
ing is correlated with fragility. This leads to 
a cross-over of the /3 values for different S at 
intermediate T as shown in Fig [Ijc). The /? 
values in Fig [IJc) are obtained by KWW fit 
to Fs{kmax,t)- However, these cross-overs are 
independent of k values as shown in Fig [T]Jd) . 
The interplay between the T-independent in- 
trinsic heterogeneity (due to the particle size 
and mass distribution) and the dynamic het- 
erogeneity which builds up at low T seems to 
be the microscopic origin of the anti-correlation 
between fragility and stretching at high T and 
the observed crossover at intermediate T. 



Fragility is usually correlated to the stretch 
exponent (3 which is found to be valid for many 
materials [T^. From PEL perspective, fragile 
liquids display a proliferation of well-separated 
basins which result in a broad spectrum of re- 
laxation times leading to stretched exponential 
dynamics [lO|- The correlation is also consis- 
tent within the framework of coupling model 
(CM) [13 according to which the strength of 
the intermolecular coupling is given by (1 — /3). 
The rate of growth of intermolecular coupling 
with decrease in T is a measure of fragility 
which according to CM would depend on the 
rate of fall of (3 with T. We indeed find that 
as S increases (fragility decreases) the rate of 
fall of /3 with T decreases (Fig[ljc)). However, 



To investigate this point in further details, 
we study the non-Gaussian parameter, a2{t) 
which also shows a correlation with fragility 
for most materials |22l |. The non-zero values of 
a2 {t) in a monodisperse system is purely due to 
the presence of dynamic heterogeneity whereas 
in polydisperse system, in addition to dynamic 
heterogeneity, there is an intrinsic heterogene- 
ity due to particle size and mass distribution 
which is present at all T. Thus for the lat- 
ter, a2{t) reflects a coupled effect of both these 
heterogeneities. As seen in Fig SI for a polydis- 
perse system a2 {t) is nonzero both in the short 
time limit (due to the mass distribution [20[) 
and in the long time limit (due to the spread 
in diffusion coefficients with particle size and 



mass). At high T, the non-zero value of a2{t) 
is predominantly due to the intrinsic hetero- 
geneity and thus increases with S (Fig^d)). 
As T is lowered, the effects of dynamic hetero- 
geneity starts to dominate, as was shown by 
the onset of connected clusters of fast moving 
particles [j, [2l| whose size increases as one ap- 
proaches glass transition. Since the relaxation 
time increases with decrease of S (Fig. [2jb)), 
there is a faster build-up of dynamic hetero- 
geneity at lower S which leads to the observed 
crossovers (Fig|31[c)&(b)) in the values of a2(<) 
between different S (similar to that observed 
for (3 in Fig[ljc)). Hence at low T, one gets the 
scenario where a2{t) decreases with polydisper- 
sity (FiglD^d)). Since fragility decreases with S, 
these crossovers in (3 and a2 (t) would mean that 
fragility is correlated only to the dynamic het- 
erogeneity and not to the intrinsic heterogeneity 
in the system. 

When the polydispersity is increased at con- 
stant volume, we get results that are opposite 
to that obtained from constant volume fraction 
studies. We find that the dynamics slows down 



with increase in polydispersity (20| . This is be- 
cause at constant volume as polydispersity in- 
creases, the packing fraction increases [231 and 
hence we find a coupled effect of polydispersity 
and density. 

Our results show that at constant volume 
fraction, although the increase of polydisper- 
sity leads to a more efficient packing, the dy- 
namics become faster due to the lubrication 
effect. Fragility decreases with polydispersity 
and is found to be correlated only to the rate of 
growth of dynamic heterogeneity and not to the 
intrinsic molecular heterogeneity in the system. 
These results reveal that the rich dynamics of 
the polydisperse system can lead to new relax- 
ations mechanisms that deserve further study. 
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